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STRUCTURAL  OF  FLAME  -  Op  UNDER  PRESSURE 


FINAL  REPORT 


The  contract  between  the  ERO  London  and  Padua  University  -  Institute 
of  Chemical  Plants  refers  to  "Studies  of  the  structure  and  composition  of 
a  flat  hydrogen-oxygen  flame  under  pressure  by  means  of  thermocouples  and 
spectroscopic  methods.  The  experimental  results  will  be  mathematically 
correlated" • 


ABSTRACT 


The  report  illustrates  the  work  carried  out  over  5  years,  suuiraarizing 
and  correlating  tb^t  which  has  already  been  set  out  in  eight  technical 
reports  on  the  state  of  progress  of  the  research  with  the  addition  of  how 
much  has  been  done  since. 

Particularly  described  are:  the  reactor  and  subsidiary  equipment, the 
equipment  and  measuring  methods  of  temperature  profiles,  concentration  of 
radical  species,  combustion  velocity  and  finally  a  new  hybrid-stochastic 
method  of  flame  simulation. 

Particular  emphasis  is  given  to  how  much  has  newly  been  done  from  an 
experimental  and  theoretical  point  of  view  and  to  the  results  obtained  in 
a  field  "flame  under  pressure",  heretofore  s tudied  1  ittle , and  in  v;hich  the 
results  given  in  literature  are  often  contradictory. 


INDEX 


1.  -  experimi:ntal. 


2.  -  MATUEMATICA].  TREATMENT. 


3.  -  rJ-.SlJLTS. 


CONCldJ.SiONS. 


RKl'i'tKRtl-r. 


—  THCRV.OCOUPlE 


THERMXOUPLE 


I 


CAPItLARY 


Fig*  4  -  Gas  flow  sifsicm. 


StPE  RAIOR 
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r/aCK  ra.  11  LC'iKowic  Dr.viCiiS 


EXPERIMENTAL. 
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a )  Reactor  and  subsidiary  equipment. 

The  lay-out  of  the  reactor  built  by  us  and  used  in  the  study  of  flames 
under  pressure  is  p,iven  in  Fig.l. 

The  characteristic  parts  are: 

-  the  burner  with  corrugated  matrix;  which  ensures  a  perfect  mixing  of  the 
entering  mixture  and  laminar  plug  flow  in  the  reactant  mixture. 

-  the  combustion  chamber,  ^  double  transparent  quartz  tube 

that  allows  working  pressure  uf  up  to  6-7  10^  Pa  and  the  carrying  out  of 
optical  and  spectroscopic  measurements  and  observations  of  flame  geometry 
and  on  the  position  of  probes  and  thermocouples. 

-  the  flameholder  is  a  holed  double  disk  in  steel  and  pyrophyllite. 

The  combustion  cliamber  is  equipped  with  automatic  scanning,  Fig. 2  wich 
allows  on  one  hand  a  high  spatial  resolution  in  the  measurements  by  means  of 
a  thermocouple  introduced  in  the  flame  as  well  as  in  the  optical  and  spectro 
scopic  ones, and  on  the  other  velocity  in  the  experimental  surveys  thus 
ensuring  constancy  in  the  behaviour  of  the  flame. 

Fig.  3  shows  the  arrangement  of  the  burner  and  instruments  both  for 
combustion  speed,  flame  temperature  and  concentration  of  chemical  species. 

Fig.  4  and  5,  in  particular,  show  in  detail  the  gas  flow  system  and 
the  optical  instrumentation . 

The  various  pieces  of  equipment  are  described  in  the  paragraphs  that 
concern  the  control  of  the  progress  of  the  reaction  or  measurement  particulars. 

b)  1'  ]  o  v;  Velocity  and  Pres  5^  u  r  e  Me  a  s  u  r  o  m c  n  t  s  , 

Next  to  the  reactor  a  critical  orifice*  is  set  up  which  produces  a 
prcssuri*  drop,  so  that  the  pressure  in  the  combustion  chambcT  is  maintained 
at  a  ctnistant  vaUie.  Critic/tl  oriiicc*s  are  also  placed  in  llu*  feeding,  lines, 
these  allow  the  flow  velocity  of  the  single  components  and  the  coropositien 
of  tin?  mixture*  to  t))e  burner  be  me  a. red.  'i’he  ])i'es5un'e  before  and  after 
is  measured  by  transducers  of  ibe  s  tra  i  n-*gaug,e  type  (Philips  Type*  PK 
set  up  ill  Whe.i !  s  I  > 'in*  bridg.i's  t>n  a  thim  sheet  of  sit'd. 

Tlie  sig,nal:;  of  tbe  lliii*e  feed  iug,  lines  and  thal  at.  t.lw  m::s  r  mil  lei 
are  me.isured  on  a  maggietir  "'i’spt'*'  s  i  i'lu )  I  .nieous  1  y  wltli  the  I  ei))}*-.  ra  i  ina*  and 
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pressure  values  in  the  combustion  chamber.  The  signals  are  also  visualised 
on  milliameters  set  up  on  the  operative  bench. 

The  transducers  can  visualise  variations  of  few  ram  of  Hg  in  normal 
working  conditions  up  to  10.10^  Pa. 

c )  Combustion  velocity. 

It  is  known  how  this  parameter  is  one  of  the  most  important  in  flame 
study.  Data  in  the  literature  are  scarce,  in  particular  for  combustion  under 
pressure.  Experiments,  if  conducted  with  accuracy,  are  the  most  trustworthy 
to  obtain  data  on  the  combustion  velocity  of  the  mixture  under  examination. 

The  measurements  are  conducted  visualising  the  flow  by  means  of 
suspended  particles  of  MgO.  The  stroboscopically  illuminated  particles  arc 
photographed  every  two  seconds  with  0.5  seconds  exposure  time.  The  camera  is 
an  Asahi-Pentax  with  a  135  mm  F  2.5  telephoto  lens.  A  double  enlargement  is 
obtained  and  a  focal  range  from  1  to  2  mm.  The  range  is  focused  on  the  luminous 
particles  that  are  found  on  the  axial  diameter  of  the  flame. 

The  velocity  is  determined  by  diagramming  the  height  of  the  particles 
against  time.  The  tangent  to  the  curve  at  a  point  gives  the  combustion 
velocity  at  that  point, 

d )  Temperature  Profiles. 

Tiic  temperature  profile  of  a  flame  is  that  wliich  analyses  more  conip]ctely 
its  evolution,  visualizing  the  preheating,  reaction  and  post-reaction  zones. 

The  temperature  values  along  the  flame  are  obtained  through  iaeasureiu<-nts 
by  means  of  thcrmoinctric  probes  immersed  in  the  f Janie  and  of  SjK'cial  diirenrions 
and  sliape,  in  order  to  disturb  the  gas  flow  as  little  as  possible. 

The  probe  used  is  a  25  microns  diameter  chromcl-alumc']  thoianocouol e 
covered  wi tli  a  fused  mixture  cl  borax  and  alumina  to  oJiminate  the  caLal3*Lie 
effects  as  n:u c li  as  p o s s i b  1  e ,  'i'  1  le  very  sma  11  d i ame t  e r  a n d  t.h e  a u t c»ma  t i  c 
scanning  allow  foi'  good  spatial  resolution. 

The  instant  position  of  tlu*  prolK'  is  recordotl  ])y  the  aid  of  a  jm>vin;; 
transducer  wliich  is  capable  C'f  resolving  the  positlcM*?  with  ^  10  microns 
nccurac)  .  Tlu‘  effective  nieasureiiicnL  time  i  .s  a  fe\v^  sce:onds,so  as  to  onsuvi 
the  coiUitancy  of  the  e»peraliv'a  vairiablt'S.  Owing  to  the  hrieO'  .‘Hi  i  C'mv’u  (  t  !i.  •, 
it  is  iv't  })ossil)Ic  li'  use  a  rioii'.l  pen  ricoixUn* ,  wlmse  cuns.  tanl  is  hig.lnr 
than  M)ai  ('f  I  ))>’  HJta.'urciiM'ji  |  ,  ihtoi'fcne  a  n'agjie  l  i  c  tape  irc'oint  r  (An.iJo;-.  7 
Philij)'.)  is  inu’d  ,  Tin-  nil  live  Ms  ar*.-  t  ra  :ai  I  im  i:n'd  ilinnipji  a  is>  1  ynan.i  .i  1  id 
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the  7th  order  which  fits  the  data  of  the  readings  and  the  temperature 
obtained  is  then  corrected  for  radiation  losses  of  the  thermocouple  junction. 

e )  Concentration  measurements  of  the  OH  and  H  species. 

The  concentration  values  of  the  radical  species  are  an  indispensable 
complement  to  the  temperature  measurement  and  a  useful  indication  on  flame 
behaviour,  in  particular  in  the  reaction  zone  where  the  pressure  effects  can 
have  a  great  effect  on  the  reaction  kinetics. 

The  relatively  lovj  flame  temperatures  allov/  the  appli¬ 

cation  of  the  emission  spectroscopy  except  in  tlie  case  of  the  fundamental 
lines  of  alkaline  metals  injected  in  the  flame.  This  technique  permits  the 
determination  of  the  H  concentrations.  The  emission  intensity  of  the  sodium 
D  lines  is  due  only  to  the  chemiluminescence  according  to  the  reactions: 

H  +  H  +  Na  - ^  +  Na 

H  +  OH  +  Na - P-  H^O  +  Na 

since  the  thermal  emission  is  irrelevant,  because  of  the  low  temperatures 
of  the  flame . 

On  the  other  hand,  absorption  spectroscopy  is  used  for  OH  concentration 
measurements.  For  such  measurements  a  quartz  xenon  lamp  is  used  as  light 
source  and  a  Philips  1003  photo-multipl  ier  as  transducer. 

To  raise  the  spatial  resolution,  cylindrical  lenses  focusing  on  the 
centre  of  the  flat  flame  arc  used.  The  image  on  the  entrance  of  the  spectrometer 
a  (Jarrel-Ash)  is  enlarged  tv;ice.  The  Oil  spectrum  band  (0,0)  from  3060  to 
3200  A,  is  taken  in  the  2nd  order  to  raise  the  I'e solving  power  of  the 
spcctromc  ter. 
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2.  -  MATHEMATICAL  TREATMENT. 


a)  Heat  Flow. 


The  volumetric  chemical  heat  release  rate  results  at  any  point 


d  dt 

—  (ra  c  T  -  X  — ) 
dz  p  dz 


(1) 


m  =  p  V  ~  mass  flow  rate  of  the  mixture  (p  ,  v  density  and  burning 
^o  o  o  o  ® 

velocity  at  normal  conditions) 

Cp,  X  =  specific  heat  and  local  thermal  conductivity. 

Equations  (1)  describes  the  heat  release  rate  as  the  gradient  of 
energy  flow  due  to  the  mass  transfer  (convection)  and  to  the  heat  conduction. 
It  does  not  make  into  account  the  heat  losses  which  are  due,  in  particular, 
to  the  radiation.  The  equation  has  been  numerically  calculated  by  Stirling's 
method. 


The  integrated  heat  rate  througouht  the  flame  interval 

z  = 


%  = 


q  dz 


(2) 


z  =  0 


can  be  compared  (apart  from  the  small  heat  losses)  with  the  heat  calculated 
from  the  reaction  enthalpy  A 

%  •  P'o  ''o 

V  >  V  =  volumetric  flow  rate  of  and  total, 
o  tot  2 

A  comparison  between  the  two  heats  allows  good  agreement  with  Ihe 

temperature  measurements. 


b)  Stochastic  h y b r i d  model  of  flame. 


In  the  mathematical  simulation  of  llio.  burning  process  of  a 
monodimens iona]  flame  stochastic  solution  lias  bec'n  adopted  ftn'  llu;  clu'inicaJ 
procofis,  which  has,  a  more  consistent  physical  ba;  is  than  the  commnr.ly 
used  deterministic.  It  adapts  better  lo  tli(‘  use  of  tin*  computer  (avo  i  lul  i  ir- 
"round -oi  l. *'  errors)  and  tlir'  ust'  of  ealculating  nKsans  to  ri  ndei*  tlie 

ado{)tec!  model  more  c*xact  and  eoncira*. 

Tlie  stc)('asli('  soJnl.iou  jn*esents  a  guod  elaslieity  sc'  miuh  5;o  Mi;  1  it 
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has  been  possible  to  hybridise  it  with  deterministic  solutions  for  energy 
transfer  and  molecular  diffusion. 


Volume  V,  spatially  homogeneous  and  in  thermal  equilibrium,  contains 
a  mixture  of  various  chemical  species  which  react  among  themselves  in 
different  ways.  They  are: 


N.  =  number  of  particles  of  the  various  chemical  species  S^: 


1  “  1,2, ••••«••, «*,n. 

=  chemical  reactions,  r  =  1>2, . ,ra. 


The  probability  of  the  reaction  taking  place  is: 

=  Qj,/  i  Q. 

1 

and  the  probable  time  is: 


(4) 


ra  m 

P„(t)  =  Q_  exp  -  (2.  Q  t)  (5) 

1  r  ^  Tc 


where  Q  =  c  f  ,  with  f  combinatory  function  of  N.  and  c  ,  stochastic 
r  r  r  r  -  i  r 

•  a— 1 

reaction  constant  =k^q!N^  ,q=  1,2,3  for  differents,  two  or  three 

equal  reacting  molecules  and  g  =  reaction  order. 

In  order  to  select  the  event  and  the  probable  time  it  is  necessary 
to  generate  two  casual  numbers  and  included  between  0  and  1  and  to 
choose  r  so  that: 

r-l  ra  j: 

1  1  1  ^ 


and 


t  =  (U  L  Q)  In  (1/L  ) 
1  ^ 


(7) 


Once  r  is  known  the  values  of  are  updated  at  time  t  according  to  the 
stoichiometry  of  reaction 

Once  the  reaction  enthalpies  and  the  heat  capacities  of  the  varioics 
componcp.ts  are  kntn^^n,  we  have  for  every  calculation  cycle  s,  and  for  the 
adiabatic  flame: 

T  =  T  -  (All)  /Z  (c-  N./N)  .s  =  0,1,2,'.... 

S  S  J.  j.  .  pi 


and  (hi‘  av«M'ap,e  fJ(>w  veloiMly  will  be; 

(NT) 


v  ^  V 


s-1  (NO 


0 
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The  conservation  equation  for  the  species  i  in  the  case  of  "premixed, 
stationary,  laminar,  ideally  nionodimensional  flame"  is : 


d 

dz 


(v-^Vd) 


■] 


dn* 
_ 1 

dr 


(8) 


with  z  =  axial  coordinate,  n  =  mole  density  number,  v  and 
=  mean  velocity  and  diffusion. 


Introducing  the  fractional  mass  flux  we  have: 


P  V 
o  o 


M. 

1 


d  G. 
1 


=  d  n. 
1 


(9) 


with  p  V  -  mass  flux  velocity 
o  o  .  ^ 

=  weight  fraction  of  species  i  in  the  flux, 

M  =5  molecular  weight. 

Since  the  variation  of  takes  place  in  a  stochastic  way,  per 
unitary  increase,  the  displacement  of  one  molecule  will  happen  when: 


p  V  G,  Z'  >1  negative  diffusion 

-2_£  _i  I  (10) 

V  1.  <  -1  positive  diffusion 

For  the  type  of  flame  adopted  the  energy  conservation  equation, 
integrated  with  the  boundary  conditions  z  =  0,  T  =  T^,  3T/9z  =  0  offers 


AT 

Az 


V  p  c  {T  T  ) 
o  o  p  o 


(11) 


with  and  X  -  mean  specific  heat  and  average  conductivity  of  the  mixture. 

Equation  (11)  is  called  upon  only  for  unitary  events  linked  with  the 
initial  reactions  between  molecules  (reactions  9  and  10  of  Table  2).  When, 
these  reactions  which  arc  the  most  probable  in  the  high  part  of  the  flame 
arc  defined  by  the  programme,  their  stochavStic  time  determines  the  entity 
of  heat  transfer. 

Ill  this  way  both  for  molecular  diffusion  as  well  as  lieat  transfer  a 
system  of  hybridization  is  chosen  adopting  llu!  deterministic  ccjuatiuns 
only  for  unitary  events  linked  to  the  stochastic  approach. 

Thus  it  is  possible  to  calculate': 

-  the  nninber  of  events; 

--  the  le.mporature  aiul  co\u‘enM-.\t  i('n  pvoliles; 
tlu;  s  I  lU'lia:;  I  i  e  t.ii'ie  veloeily  aiirl  liie  dirlancv'. 
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RESULTS. 


Among  the  numerous  runs  described  in  the  8  reports,  two  have  been 
chosen  which  characterize  better  the  various  theoretical  and  experimental 
aspects  and  on  the  other  hand  are  sufficient  in  describing  the  behaviour 
of  a  flat  flame  under  pressure. 

The  characteristics  are  given  in  Table  1. 


TABLE  1 


Chamber 

pressure 

Ingress 

Composition 

O2 

Flow 

Velocity 

Qu " 
eq(3) 

Qs- 

eq  (2) 

% 

•T. 

cm  sec 

1 

16.15  5.40  78.45 

9.05 

2.80 

2.55 

3 

12.40  5.15  82.55 

12.85 

11.30 

10.51 

, Characteristics  of  flame  at  atmospheric  pressure  and  under  pressure. 

As  regards  the  heat  release  in  general,  a  discrepancy  of  almost  10% 
is  found  between  the  values  calculated  by  equation  (3)  and  those  deduced 
by  the  temperature  profile  equation  (2).  The  agreement  is  good,  taking 
into  account  the  the  heat  due  to  radiation  is  not  calculated. 

The  reaction  probability,  the  temperature  profiles  and  the  composition 
of  the  two  flames  have  been  calculated  with  the  help  of  a  Laben  70 
computer . 

It  is  assumed  N  =  50.000. 
tot 

Of  the  38  reactions  that  can  happen  in  flames 

necessary  to  obtain  satisfactory  results  from  tlu'"  simulation  without  making 

the  calculation  too  invtU ved .The  10  reactions  are  given  in  Table  2  togetlun* 

with  the  kinetic  constants.  These  last  are  suggested  by  DAUI.CIl-DUVSDALE  - 
2 

!I0^^^•>LL0YD  (  ). 

Reaction  enthalpies,  specific  befits  and  thermal  conductivities  arc? 

3 

those  given  in  the  tabic  by  rkKNldl  (  )  ('xcej^t  ttiat  foi*  the  values 

given  in  ref.  (  )  have  hcvn  used. 
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TABLE  2 


Reaction 

,  .  ,3  ,-l  -1, 

k  cm.  (cm  mol  sec  ) 

AT,K 

AH  298(kJ/mole) 

1) 

O^+H 

O+OH 

2.2*10^^  exp(-8A50/T) 

700t2200 

70.672 

2) 

H2+O 

H+OH 

1.8«10^°  T  exp(-AA80/T) 

A 00 T 2000 

8.255 

3) 

U2+OH 

E^O+H 

2.2*10^^  exp(-2590/T) 

300t2500 

-63.30A 

A) 

H+HO2 

OH+OH 

2.5*10^^  exp(950/T) 

290v800 

-159.979 

5) 

H2O+H 

H^+OH 

9.3*10^^  exp(-10250/T) 

300-:2500 

63.304 

6) 

n+H+N, 

r  "2*"2 

2.5*lo’-®/T 

300 f5 300 

-A35.973 

7) 

O^+tt+N^-^HO^+N.^ 

1.5*10^^  exp(500/T) 

300t2000 

-197.096 

8) 

H+HO^ 

H^O+O 

5*10^^  T°'^  exp(-2000/T) 

300t350 

-231.538 

9) 

«2"02 

^  OH+OH 

S'lO^'^  exp(-22500/T) 

298t2000 

78.927 

10) 

H2^02 

H+HO2 

5.5*10^^  exp(-29100/T) 

298v2000 

238.906 

The  diffusion  velocity,  owing  to  the  large  excess  of  in  the 
flame,  is  approximated  by: 

ND.  „„  d  N. 

(v  ).  - - LiH _ 1 

'•  d'i  N.  dz  N  ' 

J  ,  - 


The  binary  diffusion  coefficient  is  evaluated  from  the  well-known 


equation. 


D.  „  =  2.628 
1.N2 


7  % 

oT.  P  Q  (T  ) 
^3 


—  10  ^  cm*'  sec  ^  (13) 


The  values  in  the  constants  are  given  in  ref.  (  )  while  the 
corrections  caused  by  the  presence  of  polar  species  arc  valued  according 

to  (  )•  In  the  simulation  programme  the  reduced  collisional  integral  H 

*  .  .  .  .  .  “1 /2 

(T  )  is  fitted  with  a  10th  degree  polynomial  expression  in  (T)  from 

7  .  .  .  .  . 

the  data  in  (  ) ; the  polynomial  coefficients  are  given  in  Table  3 

TAnidC  3 

b  1  .036917  c  =  “0.036973 

d  =  0.818603  e  -  -1.416093  f  -  2.133107 

g  =  6.38690]  h  -  -2.910819  1  -  -7.933206 

1  -  -3.843663  m  9. 390/01 

The  prog.r.ii.iikv,  are  ;;Lven  in  the  ap)>(.'rid  i  \ . 
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a)  Contribution  of  every  single  model  reaction. 

The  calculation  programme  takes  into  account  the  stochastic  path  in 
the  form  of  number  of  event  of  the  various  elementary  reactions  carried 
out  by  the  same  programme  (Fig,6),  Such  data  have  allowed  the  programme 
optimization  to  reduce  notably  the  troublesome  kinetic  model  discaring 
the  radical  reactions  which  in  our  case  were  practically  nil. 

From  the  diagrams  of  Fig. 6  it  is  possible  to  observe: 

-  the  reactions  between  radicals  tend  to  diminish  at  the  approach  of 
equilibrium.  This  can  be  represented  by  reactions  (3)  and  (5) . 

3 

+  OH  H^O  +  H 

5 

that  reach  in  proximity  of  equilibrium,  intervention  values  near  to  100%. 

-  reaction  (9)  between  H2  and  O2  tends  towards  zero  corresponding  with 
the  ignitions  point  of  the  flame. 

“  reaction  (4)  shows  a  maximum  at  intermediate  temperatures,  a  zone  in  which 
the  species  H  O2  is  stable,  and  which  later  decreases  to  zero  at 
equilibrium.  ' 

^ )  Temperature  Profiles. 

The  temperature  profiles  are  given  in  Fig. 7.  A  dotted  line  shows 
the  experimental  ones,  a  broken  line  indicates  the  calculated  ones.  The 
agreement  can  be  regarded  as  satisfactory  taking  into  account  that: 

-  the  experimental  values  suffer  from  approximate  measurements  and  an 
imperfect  flat  shape  of  the  flame. 

-  in  the  model  the  flame  is  assumed  adiabatic  and  certain  secondary  reactions 
are  disgrcgarded. 

From  the  two  profiles  of  1  and  3  atm  it  is  possible  to  note,  in 
agreement  with  the  increase  of  the  global  reaction  velocity  with  j) re a sure: 

In  (rj^/r2)  =  j^(q/2)  ij  In  (p^/p2)  v’hcre  q  is  tl)C  global  reaction  order  - 
=  2,081  (^),  that  equilibrium  is  reached  more  q\iickly  as  tint  pressure  is 
increased  and  the  reaction  zone  becomes  thinner. 

It  is  als{'>  noted  (ruu5i  des(‘ribed  in  p)*evj(ms  reports),  tliat  the 
effect  of  diffur.ion  is  inorc’  accc'ut  u.i  (  tcl  at  lev;  pri'j.rairt'  than  at  high 
prer.sina;,  givcui  that::  tla*  g,lo]>al  react  ion  V'‘lc>cjly  is  proport  ituril  lo 
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p(^/^)“l  ;  the  diffusion  is  inversely  proportional  to  pressure,  as  it 
increases,  the  reactions  gain  in  importance. 

c)  Concentration  Profiles. 

The  facts  previously  described  find  verification  also  in  the 
concentration  profiles,  Fig.  8-9.  The  increase  of  the  global  velocity 
leads  to  a  quicker  reaching  of  equilibrium  conditions. 

The  two  most  important  chain  carriers,  OH  and  H,  show,  as  intermediate 
species,  concentration  maxima  whose  position  depends  on  pressure.  The 
relative  curve  at  OH  presents  two  of  them,  the  first  relative  to  reaction 
(4)  masimum  and  the  second  together  with  that  of  H,  due  to  reactions  (1) 
and  (2)  maximum. 

The  presence  of  a  minimum  for  OH  had  been  noted  in  the  experimental 
determination  of  such  species,  report  6. 


CONCLUSIONS  . 


The  theoretical-experimental  study  of  flames  pressure 

has  allowed  us  to  identify  the  influence  of  the  pressure  parameter  on 
fluid  dynamics  and  on  the  progress  of  the  system  reactions  in  flow. 

The  small  reaction  zone  of  flame  under  pressure  and  the  low  temperature 
of  flame  due  to  the  excess  of  are  favourable  conditions  in  the  assumption 
of  single  dimensional  flow  and  therefore  ideal  for  the  experimental  chech  of 
the  simulation  model  adopted. 

This, a  hybrid-stochastic  type,  offers  better  oppor tuni ties  in  respect 
of  the  deterministic  normally  used.  In  fact  besides  avoinding  "round  off 
errors",  giving  stable  and  accurate  solutions,  it. is  capable  of  indicating 
the  most  important  physical  phenomena  and  the  most  probable  reaction 
mechanism  for  the  various  temperature  and  pressure  conditions  of  the 
flame  considered. 
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